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Abstract 

We consider the problem of pricing convertible bonds deriving an analytic solution for the value of 
an option to convert the bond to equity at the bond’s face value at bond maturity. The formula derived 
is based on the perturbation expansion techniques developed by Turfus (2016a) and Turfus (2016b). 
The main result is an analytic expression for the correction to the value of the option arising from 
the expected negative correlation between the equity price and the credit spread for the bond issuer. 
Somewhat counter-intuitively the correlation risk is seen to be wrong-way in the sense that it decreases 
the value of the embedded optionality. Although the method assumes that the impact of the credit 
intensity is weak, it is proposed that it ought to give reasonable approximations even for realistic market 
values. 


1 Introduction 

We look to price the embedded European option associated with a convertible bond which allows the holder 
to exchange the bond for a fixed number of the issuer’s shares at maturity, in lieu of repayment of notional. 
Our model is closely related to one proposed by Davis and Lischka (1999), specifically the variant in their 
§2.3 where interest rates are taken to be deterministic but hazard rates stochastic. They consider the more 
general problem of Bermudan options to convert during the life of the bond and adopt a tree-based approach 
to facilitate the calculation. They use a Hull and White (1990) short rate model to enable greater tractability 
in their numerical computations, whereas we prefer to adopt a Black and Karasinski (1991) model, ensuring 
as it does that hazard rates remain positive. While they take into account correlation between the equity 
and the default process, they do not take into account the default (jump) risk on the equity process which 
results in an increased drift, as we shall see, when the option to convert is priced contingent on no default 
having occurred. Later authors such as Ayache et al. (2003) who look to handle the default contingence of 
options to convert more carefully similarly neglect the modification which is needed to the equity dynamics 
when it is modelled contingent on no default. 

We set out in section 2 below a description of the underlying stochastic model for the equity and default 
intensity processes, going on to show how a PDE describing the value of an option to convert a bond to 
equity at face value can be derived. In section 3 we show how, under an asymptotic assumption of a weak 
default intensity, an analytic expression for this can be inferred as a perturbation expansion. Our main result 
is expressed in the concluding section 4, in particular in Theorem 4.1. 


2 Model Description 

We consider a credit default process driven by a stochastic default intensity, the dynamics of which are 
governed by a Black-Karasinski short rate model. The equity price we take to be given by a jump-diffusion 
process, with a downward jump to zero at the time of default. The diffusive processes are mutually correlated. 
We follow Turfus and Schubert (2017) in their related calculation for contingent convertible (CoCo) prices 
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in assuming from the outset that rates are deterministic. Thus we propose: 


^ = -a\(t)(\nX t -\n6\(t))dt + ax(t)dWl, (1) 

A* 

= ( r(t ) - q(t) + X t ) dt + a s (t ) dW? - dn t . (2) 

for t < T m = min(r, T), where r is the time of default of the bond/equity issuer and T is the maturity of 
the convertible bond under consideration. Here A t is the default intensity with mean reversion rate a\(t) 
and mean reversion level 9\(t) chosen to satisfy the no-arbitrage condition that 


E 


e - fo ds 


= e - fo A ( s ) ds 


(3) 


under the martingale (money market) measure for 0 < t < T, with A (t) the instantaneous forward default 
intensity observed at the initial time t = 0. Here St is the equity price, r(f) is the instantaneous forward rate 
of interest, q(t) is the expected (continuous) dividend rate, dW} and dW/ are correlated Brownian motions 
with 

corr(W t \ W t 2 ) = p XS 

and rit is a Cox process with intensity A t giving rise to an equity price jump to zero, contingent on default. 

We suppose that the convertible bond is a standard bond priceable by the usual means, with an embedded 
option to exchange the notional repayment at T for a fixed number M of ordinary shares per unit notional. 
We write the PV of a convertible bond with notional N as 


V{t) = Nf hond (t) + Nf option (t). (4) 

where /bond (i) expresses the part of the PV generated by bond cash flows, i.e. the coupon payments 
and notional repayment, and /option(t) expresses the additional contribution to the PV resulting from the 
embedded option. We focus here on the calculation of the latter term. Making explicit the dependency 
on the market variable S t and A t , we denote /option (t) = We deduce by standard means that 

f(S, A ,t) will be governed by the following backward diffusion equation: 

§[ + + ^ + 2 P\s<^s<J\SX-^^ + a 2 x X 2 ^^j 

= 0 f(t ) + A)/ 

(5) 


for t < T m and subject to the assumed terminal condition that 

/(St, X t, T) = rri ax {A# SV — 1,0} 

with M the number of shares issued per unit notional. We seek to calculate /(<So,Ao,0), the value at time 
t = 0 of the embedded option to convert. 

3 Perturbation Analysis 

In the absence of an exact analytic solution to (5), we seek an approximate solution. The perturbation 
approach used here to obtain convertible bond prices was first proposed by Turfus (2016a), who proposed 
that the short rate in a Black and Karasinski (1991) model be considered weak in an asymptotic sense. It 
is closely related to that used by Turfus (2016b) in calculating the prices of CoCo bonds under an assumed 
Black-Karasinski short rate model for the conversion intensity, combined with an equity jump-diffusion 
process. This was an extension of the earlier calculation of Turfus and Schubert (2017) based on an assumed 
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Hull and White (1990) short rate model. However in both these CoCo bond calculations an assumption was 
made that the credit intensity volatility was small, rather than the credit intensity itself, as proposed below. 
Particularly in the Black-Karasinski case the former approximation leads to less accurate results, whence we 
suggest the approach set out below is superior. 

Following Turfus (2016a), we suppose that the default intensity is weak, i.e. that an asymptotic parameter 

e = fpXt)dt 

Ip a\{t)dt 


is small, wherefore we define a scaled forward intensity A(t) = ^A(f). We proceed by noting that the solution 
to (1) can be written 

In At = lnA*(f) + f (j>(s,t)cr\(s)dW 
Jo 

where 

:=e-£ a ^ u)du (6) 


for some as yet unknown O(e) function X*(t). For computational convenience we propose a new stochastic 
intensity coordinate 


y t := In A t - In A*(f) + 


(7) 


with 


We also write 


I\(ti,t 2 ) := [ <p 2 (v,t 2 )<J 2 \(v)dv, 

Jt i 


X*(t)=eXl(t) + e 2 X* 2 (t) + 0(e 3 ) 


( 8 ) 


where we know from Turfus (2016a) that the no-arbitrage condition (3) and the initial condition yp = 0 
require us to choose 1 


a m 

= A(i) 

(9) 

m 

= A(t) [ (. F x (v , t) - 1) X(v) dv. 

(10) 


o 


with 

F\(v,t) := exp t)I\(0, v)). (11) 

Similarly, we define a new characteristic stochastic equity process Xt such that 

S t = F(t)£(x t ), (12) 

where £(X t ) := exp (X t - \[X] t ) is a Doleans-Dade exponential with [X\ t the quadratic variation, and 

F(t) := Sgeti^-q^+X^dS' ( 13 ) 

We define also for notational convenience an equity payoff function 

M 0 (x,t) = MF{t)e x ~i Is( - °’ t \ 

where 2 t 

Is(ti,h)= [ cr%{u)du. (14) 

Jt i 

1 We will not in fact use the expression for to the level of approximation explored in the present context. 

2 Note that ls(0,t) is the quadratic variation of xt required for (12). 
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Changing to the new coordinates (x, y , t) and expressing 


(5) can be rewritten as: 


with final condition 


h(x, y, t) = 


C[h] ~(A*(/)e ?/ - eA(t)) (h - §|) 


h(x T ,yT,T ) = max{M 0 (x T ,T) - 1,0}, 
here £[.] is a standard forced diffusion operator given by 

A 8 1 / r) 2 r) 2 h r) 2 \ 

= S - + 5 (^> w + + <4(0 - «0 + W). (16) 

To solve this we pose an asymptotic expansion 


At zeroth order we must solve: 


with final condition 


h{x, y, t) = h 0 (x, y, t ) + e/ii(x, y, t) + 0(e 2 ). 


C[h 0 ] = 0 


ho{x T ,y T ,T) = max{ M 0 (x t ,T) - 1,0}. 
This can be achieved by means of a Green’s function 


G(x, y, t; £, rj, v) = B(t, v)H(v - t) AT 2 (£ -x,y- ycj>(t, v);R{t , u)). 


where H(.) is the Heaviside step function, 


= e -/ t ”Ms)+A( S ))d S; 


is a risky discount factor, and N 2 (x, y, R{t, v)) is a bivariate Gaussian probability distribution function with 
mean 0 and covariance matrix 


where 


Making use of (18) we obtain 


*“•■>-( $3 $ 3 ) 

rt2 

lR{tl,h) = P\S / (/>i v R2)cr\(v)(Js(v)dv. 


where 


/ OO nOO 

/ G(x u y t , t; £, 77, T) max{M 0 (£, T), 1} d£ dy 

-oo J — OO 

= B(t, T) (£(x t )MF(T)N(d 1 (x t ,t,T)) - N(d 2 {x t ,t,T))) 

1 f x 

N(x) = . _ / exp (— |w 2 ) du, 

V 27T J —oo 


d 2 (x,t,v) = 


lnMo(x, x) — In AT 


di(x,t,v) = d 2 (x, t,v) + \/ls(t,v). 
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We note that there is no dependence on the stochastic intensity variable y t at this level of approximation. 
Continuing to first order, we find we must solve 

CM = A (t) 0,t) _ ^ ( /l0 _ , 


Proceeding as before we obtain 


hi(x t ,y t ,t) = B(t, T)£ (x t ) 



M u ) (c (: y t <j)(t,u )) 


g In{t,u) 



N(d 1 (x t + I R {t,u),t,T))du. 


(26) 


4 Conclusion 

Substituting the above expressions into (17), setting t = 0 and Xq = yo = 0 and resorting to unsealed 
variables, we obtain our main result: 

Theorem 4.1 The value of an embedded, option for the holder to convert a bond to M shares of the issuer’s 
stock per unit unit notional at maturity T is, under the assumption that the evolution of the credit default 
intensity and the share price are given respectively by (1) and (2), given asymptotically in the limit of small 
credit default intensity by 

fo P uon( 0) - R(0,T)(mF(T) (^(0,0,7)) + T -jS== (e lR ^ u) - l) JV'(di(/ fl (0,«), 0, T) du) 

V V Jo y/I S (0,u) V 7 J 

- N(d 2 (0,0,T))^j 

(27) 


with errors = where the notation is as defined in the preceding section. 

This formula is closely related to the defaultable variant of the well-known Merton (1973) formula and indeed 
reverts to it in the limiting case of deterministic credit intensity or when p\s = 0. The difference is in the 
inclusion of the integral term in (27). Somewhat counter-intuitively, the impact of the expected negative 
correlation between equity price and credit spread is seen to result in a slight decrease in the option value: 
the obvious impact of increased survival probability associated with higher-than-expected values of St is 
more than offset by the concomitant decrease in the jump-compensating drift (the X t term in (2)). 

Based on the high degree of accuracy found by Turfus (2016b) in a similar calculation using a similar 
approximation, we expect the error associated with the result in Theorem 4.1 to be at most a few bp per 
unit notional. 
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